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ABSTRACT 

Recent observations reveal that a majority of galaxies contain a central massive object (CMO), either a 
supermassive black hole (SMBH) or a compact stellar nucleus, regardless of the galaxy mass or morphological 
type. The masses of these CMOs correlate tightly with those of the host galaxies, Mcmo ~ 0.002 M ga \. Several 
recent studies argue that feedback from black holes can successfully explain the M^-ct correlation in massive 
elliptical galaxies that contain SMBHs. However, puzzles remain in spirals or dwarf spheroids that do not 
appear to have black holes but instead harbor a compact central stellar cluster. Here we use three-dimensional, 
smoothed particle hydrodynamics simulations of isolated galaxies to study the formation and evolution of 
CMOs in bulgeless disk galaxies, and simulations of merging galaxies to study the transition of the CMO- 
host mass relation from late-type bulgeless spirals to early-type ellipticals. In the simulations, absorbing sink 
particles represent either SMBHs or star clusters, while stellar feedback on the gas is implemented by assuming 
an isothermal equation of state with effective sound speed of order 10 km s" 1 . Our simulations show that the 
mass of the CMO correlates with that of the host galaxy in both isolated bulgeless spirals and in ellipticals 
formed through mergers, and that Mcmo correlates with the global star-formation efficiency in the galaxy. We 
find that the final mass of the CMO is dominated by the accreted mass, rather than the initial fragment mass. 
The denser nuclei of more massive galaxies have higher mass accretion rates, and both the final accreted CMO 
mass and the recently formed stellar mass increase monotonically with the total mass of the galaxy. Our results 
suggest that the observed correlations may be established primarily by the depletion of gas in the central region 
by accretion and star-formation, and may hold for all galaxy types. A systematic search for CMOs in the nuclei 
of bulgeless disk galaxies would offer a test of this conclusion. 

Subject headings: galaxy: nuclei — galaxy: spirals — galaxy: ellipticals — galaxy: evolution — galaxy: 
kinematics and dynamics — galaxy: ISM — galaxy: star clusters — stars: formation 



1. INTRODUCTION 

Supermassive black holes (SMBHs) appear to exist in 
most, if not all, galax ies (see, e.g., the revie ws by 
lHaiman & Ouataertll2004l and IFerrarese & Fordll2005l) . Over 
the past few years, SMBHs have been detected in about 
40 ne arby galaxies using gas and ste l lar dynamical methods 
(e.g ., Kormendv & Richstonel 11995b iRichstone et al.1 119981: 
lHoTll999t iMerritt & Ferraresdl200ll) . Thes e galaxies in- 
clude nearby ellipticals and lenticulars (e.g., iMarconi et al.1 
1 1997b Macchett o et al.lfl997l), as well as s pirals with bulges 
such as the Milky Way dGenzel et al.l 119971: iGhez et al.l 
120001 120031: ISchodel et all |2002j). Most of these galaxies 
are luminous ones with total masses larger than 10 I2 M Q . 
There appear to be tight correlations between the masses 
of the SMBHs, Mbh, and the global properties of the 
spheroid componen ts of their hosts, such as their lumi- 
nositi es and masses (Magorri an et al.lll9 98: Marc oni & Hunt! 
2003), light concentration dGraham et al.ll2001l). and velocity 
dispersions (IFerrarese & Merritlll2000F lGebhardt et al.l l2000b 
iTremaine et al.ll2002l) . 

Recently, low- and intermediate-luminosit y galaxies, with 
Mb> -18 mag, have been studied by IWehner & Harrisl 
(2006) using the Hubble Space Telescope (HST) Wide- 
Field/Planetary Camera Dwarf Elliptical Galaxy Snapshot 
Survey of the Leo Group and Virgo an d Fornax c lusters 
dLotz. Miller. & Ferguson! 120041) . and by ICote et al.l (120061) 



and IFerrarese et all (120061) . using HST Advanced Camera for 
Sur veys data for psIO Q early-type galaxies in the Virgo clus- 
ter dCote et al.l 120041) . These analyses show that a major- 
ity of these faint galaxies contain a compact stellar clus- 
ter at their center. The masses of these star clusters cor- 
relate tightly with those of their host galaxies, Mcmo ~ 
0.002M ga i, following the same relation between SMBH and 
galaxy ma sses seen in more luminous galaxies with Mr < 
-20 mag dMagorrianet all 11998b IFerrarese & Merrittl 12000b 
IMarconi & Huntll2003l) . A correlation between the mass of 
nuclear star clusters and galaxy properties holds not only 
in elliptical galaxies, but also in spira l galaxies with bulges 
dBalcells et al.ll2003t iRossa et al.ll2006l) . Nuclear star-clusters 
with masses in the range ~ 10 6 - 10 8 M ( r 1 have also been foun d 
in a handful of bulgeless spiral galaxies ( Wal cher et al.l 2005). 
These findings strongly suggest that the SMBHs in massive 
galaxies and the compact stellar nuclei in less luminous galax- 
ies may form as a result of similar physical processes, linked 
to the formation of the galaxy as a whole. 

Many models have been proposed fo r the observed 
SMBH-bulge correlations, as reviewed by Rob ertson et al.1 
(2006a), notably including self -regulation by global feed- 
back, which suppresses further accretion and star formation 
(e.g. | Silk_& Reesl 119981: lHaehnelt et ail 1 19981: iFabia^ 1 1 999b 
King! 120031: IWvithe & Loebl 12003b iDi Matteo etail l2005b 
Spring el et al.l l2005b : ISazonov et all 12005b iMurrav et al.l 
2005b Wvithe & Loeb 120051); as well as environmental reg- 
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20051) : and star formation models in which gas dissipation 



TABLE 1 

Spiral Galaxy Models and Numerical Parameters 



affect s the growth of the S MBHs (e.g., I Archibald et alj 
I2002t iDi Matteo etaT1l2003t iKazantzidis et alJ l2005h . In 
particular, self-regulated models with SM BH feedback 
of ph otoionization and Compton heating dSazonov et alj 
2005), or in the form of thermal energy coupled to the 
ambient gas have been demonstrated to successfully re- 
produce many observed properties of elliptical galaxies 
formed by majo r merg e rs, including th e Mbh- c relation 
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colors ( Sprinael, Di Matt eo. & Her nquist 
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galaxy 
X-ray gas 

emission dCox et al.ll2 006b). quasar properties and luminosity 
functions (Hop kins et al.1 2006a b), as well as the lum inous 
quasars observed at the highest redshift dLi et al.ll2006al) . 

However, self-regulated models are applicable primarily 
to galaxies with SMBHs (i.e., mass i ve ell ip tical galaxies). 
The n ew o bservations by Cote et al. (2006), iFerrarese et al.1 
(2006j) and IWehner & Harril ((2006) show that many low- 
and intermediate-mass galaxies that do not appear to have 
SMBHs still have central compact stellar clusters that obey, to 
within observational errors, an indistinguishable mass corre- 
lation with their host galaxy. This suggests that galactic black 
holes and central compact stellar nuclei, hereafter collectively 
referred to as compact massive ob j ects (CMOs), a terminol- 
ogy first introduced by ICote et all (120061) . may have formed 
from the same processes — gravitational collapse of a mas- 
sive clump of gas at the center of the galaxy. Furthermore, it 
suggests that the CMO-host mass correlation may be univer- 
sal, regardless of galaxy type or mass. 

In order to test this hypothesis, we start our investiga- 
tion of CMOs with bulgeless spirals, which are the simplest 
galaxy models, and therefore ideal to study the physical pro- 
cess of CMO formation. Nucle ar star clusters have been 
obs erved in bulge l ess spirals b y iBoker et all (120021 12004b 
and IWalcher et all (120051 120061) . These authors find that the 
growth of the CMOs is closely linked to the star formation his- 
tory of the host galaxy. We here examine the CMO-host mass 
correlation in simulations of isolated, bulgeless disk galaxies 
with a wide range of masses and gas fractions. We then further 
investigate how the CMO-host relation changes in simulations 
of mergers between spiral galaxies. 

The rest of this paper is organized as follows. In § 2, we de- 
scribe our computational methods, galaxy models and param- 
eters. In § 3, we study the formation and evolution of CMOs 
in isolated, bulgeless, spiral galaxies. We then study the for- 
mation of CMOs in mergers between disk galaxies in § 4. We 
discuss the connections between CMOs and host galaxies in 
both spirals and ellipticals, and conclude in § 5. 

2. COMPUTATIONAL METHOD 

We use the publicly available, three-dimensional, paral- 
lel, N-body /smo othed particle hydro dynamics (SPH) code 
GADGET vl.l (ISpringelet all 120011). modified to i nclude 
absorbing sink particles dBate. Bonnell. & Pricdl 19951) to di- 
rectly measure the mass of gra vitatio nally collapsi n g gas . 
iLi. Mac Low. & Klessenl (l2005bl) and Uappsen et all (2005) 
give detailed descriptions of sink particle implementation and 
interpretation. Sink particles replace gravitationally bound re- 
gions of converging flow that reach number densities > 
10 3 cm" 3 . The sink particles have a control volume with a 
fixed radius of 50 pc from which they absorb surrounding 
bound gas. They interact gravitationally and inherit the mass, 
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a Model of single disk galaxy. First number is rotational velocity in km s _1 
at the virial radius, the second number indicates sub-model. Sub-models have 
varying fractions mj of total halo mass in their disks, and given values of / g . 
Sub-models 1-3 have m L \ = 0.05, while sub-model 4 has mi = 0.1. b Virial 
radius in kpc within which the mean mass density of the halo is 200 times of 
the critical density. Virial mass of the galaxy in 10 10 MQ. d Fraction of disk 
mass in gas. c Radial disk scale length in kpc where stellar surface density drops 
by e . f Gravitational softening length of gas in pc. g Gas particle mass in units 
of 10 4 M Q . h Model run with LT mode c s = 6 km s -1 .' Model run with HT mode 
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Galaxy Mergers and Numerical Parameters 
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a Model of galaxy mergers. Note all progenitors have disk mass fraction of mi = 0.05, 
. b Rotational velocity in km s -1 at the virial radius for the first progenitor. Virial mass 
of the first progenitor in lO IO M0. d Rotational velocity in km s -1 at the virial radius for 
the second progenitor. Virial mass of the second progenitor in 10 I0 MQ. f Gas fraction of 
the progenitors. 8 Gravitational softening length of gas in pc. h Gas particle mass in units of 
10 4 Mq .' Sound speed for the gas in km s -1 . 

and linear and angular momentum of the gas they accrete. 

2.1. Galaxy Models 

Our models of isolated, bulgeless spiral galaxies con- 
sist of a dark matter halo f ollowing the prescription of 
dNavarro. Frenk. & White 119971) . and an initially exponential 
disk of stars and isothermal ga s. The galaxy structure i s 
based on the analytical work bvlMo. Mao. & White] d!998l) . 
as implemented numerically by lSpringel & W hite ( 199% and 
Springel (2000). The detailed descript ion of t he mod els and 
numerical simulations are given in ILi et all J2005bl) . Ta- 
ble 1 lists the most important model properties and numer- 
ical parameters. In order to sufficiently resolve gravita- 
tional collapse, the models are set up to satis fy three numer- 
ical criteria: the Jeans resolution criterion (Bat e & Burkerj 
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119971: IWhitworthl 1 1 9981) . the gravity-hydro balance criterion 
for gravitational softening dBate & Burkertl 1 19971). and the 
equip artition criterion for particle masses dste inmetz & White 
1997). We choose the particle number for each model such 
that they not only satisfy the criteria, but also all runs have at 
least 10 6 total particles. The gas, halo and disk particles are 
distributed with number ratio N g : Nh '■ Nd = 5 : 3 : 2. The grav- 
itational softening lengths of the halo = 0.4 kpc and disk 
hi = 0.1 kpc, while that of the gas h g varies with models. The 
minimum spatial and mass resolutions in the gas are given 
by h g and twice the kernel mass (~ 80m ? ). We adopt typical 
values for the halo concentration parameter c = 5, spin pa- 
rameter A = 0.0 5, and Hubble constant Ho = 70 km s" 1 Mpc -1 
dSpringelll2000l) . 

To test the effects of feedback, we have chosen two val- 
ues of the isothermal sound speed from the range of ve- 
locity dispersions typical ly observed for neutral gas in disk 
galax ies 6-15 km s _1 (lElmegreen & Scald I2004T iDib et all 
2006). Models designated as low-temperature (LT) have c s = 
6 km s" 1 , while high-temperature (HT) models have c s = 15 
km s" 1 . We should emphasize that the effective tempera- 
ture approach we use in our simulations should be regarded 
as a simple approximation of the underlying processes main- 
taining the interstellar turbulence level, including stellar feed- 
back, magnetorotational i nstability, and other processes (see 
iMac Low & Klessenl 120041 for a review). All of the galaxy 
models have been simulated in the HT mode. However, in 
our simulations, low mass HT models (e.g., some of the SG50 
and SGI 00 models) do not appear to collapse in the first 3 
Gyrs. Only a limited number of runs have been performed 
in LT mode, as it is prohibitively expensive to run high-mass 
LT models because they require ~ 10 7 particles to fully sat- 
isfy the three numerical criteria. The model properties and 
numerical parameters are summarized in Table 1. 

The same set of simulations of isolated galaxies has suc- 
cessfully reproduced many observations of star formation 
in disk galaxies , including distributions and morphologies 
( Li et al.l 2005bl). and both glob al and local Schmidt laws 
dLi. Mac Low. & Kle ssen 2006b), demonstrating that star for- 
mation in galaxies are controlled by gravitational instability 
(ILL Mac Low. & Kle ssen 2005a). In this paper, we focus on 
the connection between the CMO and the host galaxy. 

In order to investigate CMOs in elliptical galaxies, we 
have also done five representative models of mergers of spi- 
ral galaxies that result in elliptical galaxies. These mergers 
vary in the mass ratio of the merging pair, progenitor size, 
gas fraction, effective sound speed, bulge component, and or- 
bital parameters and configuration. Model MG4 resembles 
the merger of Milky Way and Andromeda, in which the pro- 
genitor properties and orbital p a ramet ers ar e set up fo l lowin g 
iDubinski, Mihos. & Hernquisil d 19961) and ISpringell (120001) . 
The mass ratio of the progenitors is roughly 1 :2, and each pro- 
genitor has a bulge with a mass that is 25% of the disk mass, 
but neither progenitor initially has a CMO. The two galaxies 
are on a parabolic orbit, with an initial separation of 700 kpc, 
and a pericentric distance of 5 kpc. The rest of the models are 
equal mass, head-on mergers of identical bulgeless progeni- 
tors using models SG50-1, SG100-1, SG100-3, and SG220-1, 
respectively. These mergers are run with parabolic orbits hav- 
ing pericentric distance of 1 kpc, whose separation is the sum 
of the virial radii of the two progenitors. Each of the mergers 
is run with particle number around iV to t ~ 2 x 10 6 . The model 
properties and numerical parameters of the merger runs are 
listed in Table 2. We have also performed a resolution study 



for the LT versions of models SGI 00-1 and SG220-1, running 
each with total particle numbers of 10 5 , 8 x 10 5 , and 6.4 x 10 6 , 
in addition to the standard value of 10 6 . These models cover 
a factor of four in linear resolution. 

2.2. Interpretation of Sink Particles as Clusters and Central 
Massive Objects 

The absorbing sink particles in our models trace regions 
with pressure of P/k ~ 10 7 K cm" 3 , and reach masses 
of > 10 6 M^, typical of sta r-forming proto-stellar clusters 
(Elme green & Efremovl 19971) . We therefore interpret the for- 
mation of sink particles in general as representing the for- 
mation of massive stellar clusters. We find in the simula- 
tions, however — both in the case of isolated galaxies and for 
galaxy mergers — that the most massive s ink particle alw ays 
resides in the center of the isolated disk JLT et al. 2005bl) or 
the merger remnant dLi. Mac Low. & Kl essen 2004). It either 
forms there, or migrates to the center quickly after formation 
due to dynamical friction. We therefore distinguish this sink 
particle from the others, and identify it as a CMO. 

To quantify star formation, we assume that individual sinks 
represent dense molecula r clouds that form stars at some effi- 
ciency. Observations bv lRownd & Youn3(fl9 99) suggest that 
the local star formation effi ciency (SFE) in m olecular clouds 
remains roughly constant. Kennicutt (1998J) found SFE of 
30% for starburst galaxies that lWong & Blita (120021) showed 
are dominated by molecular gas. We therefore adopt a fixed 
local SFE of 30% to convert the mass of sinks to stars, while 
making the simple approximation that the remaining 70% of 
the sink particle mass remains in gas form. This conver- 
sion rate was shown to reproduce both th e global and lo cal 
Schmidt laws observed in nearby galaxies dLi et a l. 2006b). 

We assume, on the other hand, that the CMO inherits 
the total mass of the central sink particle, for the follow- 
ing reasons. (1) If the CMO forms a compact stellar clus- 
ter, then the star formation efficiency could be > 50% due 
to extremely high pressure an d density at the galactic center 
(Elmegreen & Efremov 1997), where the deep gravitational 
potential also helps to keep the gas bound the CMO. (2) If a 
fraction of the CMO collapses into a black hole, then it could 
rapidly accrete the remaining gas. (3) Although feedback- 
driven winds from either supernovae or a quasar could expel 
the gas from galactic nuclei, the effect m ay be minor. A study 
of galactic winds by ICox et all d2006a) shows that the mass 
loss by a typical stellar wind (wind efficiency ~ 0.5, wind ve- 
locity ~ 500 km s" 1 ) is less than 10% of the initial total gas 
mass. The quasar outflow is comparable to this in galaxies 
with accreting black holes. Furthermore, most of the gas in 
the central region is consumed quickly by either star forma- 
tion or black hole accretion, so the mass loss from the CMO 
by feedback is likely negligible (ICox et al.l l2006b). 

The final mass and fate of the CMO should depend on 
the galaxy potential and dynamical evolution. The larger the 
galaxy, the more massiv e its central cl ump, due to stronger 
gravitational instability dLi et a l. 2005a). The gas clump may 
then form a supermassive black hole directly by rapid core 
collapse within a deep potential well, fragment to form a 
massive star cluster in a less-massive galaxy, or engage in 
cluster-cluster collisions to form a SMBH in ga laxy merg- 
ers due to stellar dynamical p rocesses (see, e.g., iReesll 19841 
and lHaiman & Ouataerg f2004 for reviews of different routes 
to forming a SMBH). Furthermore, the McMO _ A^gai relation 
may be established by the same physical mechanism that de- 
termines the fate of the clump — massive galaxies or galaxy 
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mergers that have higher CMO mass tend to form stars at a 
higher rate than less massive ones. These considerations sug- 
gest that the McMO-M ga \ relation is directly connected to the 
star formation efficiency of the galaxy. 

In the simulations, we do not resolve the transition within 
the sink particle from a collapsed gas clump to a black hole or 
a star cluster, and consequently we will limit our discussion, 
and comparisons to observations, to CMOs in general. 

3. CENTRAL MASSIVE OBJECTS IN BULGELESS SPIRAL 
GALAXIES 

3.1. The Formation and Evolution of CMOs 

In all of our simulations, high gas density regions collapse 
to form stars and clusters once the disk becomes gravitation- 
ally unstable. Sink particles first form in the central regions 
of galaxies; star formation subsequently propagates outwards. 
The collapsed objects continue to grow by gas accretion. 

The growth histories of CMOs in several models were 
shown in Figure [TJ The CMO initially collapses at a mass 
above the Jeans mass (M s ~ 1O 6 M and - 2 x 10 7 M Q for 
LT and HT models, respectively). We find that in the first 
few hundred million years, the CMO grows rapidly, accret- 
ing gas particles gravitationally bound to it, at a rate that is 
within an ord er of magnitude of the spherical Bondi-Hoyle 
accretion rate dHovle & Lvttletonll94UBondi & Hovlefl944l 
lBondil[T952h : 

4iraG 2 Ml MO p 



M m = 



(c 2 + v 2)3/2 



(1) 



where Mcmo is the CMO mass, p and c s are the density and 
sound speed of the gas, respectively, a is a dimensionless pa- 
rameter of order unity, and v is the velocity of the CMO rel- 
ative to the gas. After the initial accretion period, the growth 
slows down. In about 1 Gyr, the masses of most CMOs satu- 
rate due to the depletion of gas around the nucleus by star- 
formation, as show n in the middle panel of Figure [TJ We 
note that a study by ISazonov et al.l (120051) of radiative feed- 
back from accreting black holes in ellipticals shows that the 
observed Mbh-c relationship could be established once most 
(> 99%) of the gas has been converted into stars. This finding 
agrees with our results that the growth of the CMO is tied to 
the star formation in the galaxy. 

The resolution studies show that the masses of the CMOs 
converge to better than a factor of two in runs with parti- 
cle numbers > 10 6 for both models SGI 00-1 and SG220-1, 
in which just the gas particle mass is smaller than the Jeans 
mass at which the CMO collapses initially. We include in our 
analysis two models that satisfy this requirement in order to 
extend the range of LT models that we can consider, as we 
will discuss. 

3.2. The McMO—M ga \ Correlation 

In Figure [2] we show the correlation between the mass of 
CMO and that of the host for the bulgeless, isolated disk 
models. We consider three different definitions of galaxy 
mass: M sg , the combined mass of stars and gas, M tot , the total 
mass of the galaxy including stars, gas, and dark matter, and 
M ga i,disk, the total mass of the galaxy within 5R e , where R e 
is the half-light radius, containing half of the gas and stellar 
mass of the disk component. In all three cases, Mcmo is mea- 
sured after 1 Gyr, or in the few cases where the simulations 
stopped before 1 Gyr, at the last time step. The two LT mod- 
els SGI 60-4 and SD220-4 (shown in thin circles) do not fully 
satisfy the Jeans criterion (i.e., the SPH kernel mass of 32 gas 
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FIG. 1 . — Mass growth histories of the CMOs in selected galaxy models 
(top panel), and the stellar mass fraction in the nuclear region within 1 kpc 
of the galactic center (middle panel). Also shown are the histories of CMO 
mass from the resolution study (bottom panel). The designations Rl, R8, R64 
and Reg indicate total particle numbers of 10 5 , 8 X 10 5 , 6.4 X 10 6 , and 10 6 , 
respectively. 
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FIG. 2. — The relation between CMO mass and host mass in simulated bulgeless spiral galaxies. Galaxies are shown from both the LT (open symbols) and HT 
(filled symbols) models. The color and the shape of the symbol indicates the rotational veloci ty and sub-m odel (gas fraction and disk scale length) for each galaxy 
described in Table 1 . The solid line is the best fit to spheroid systems observed by Ferrarese et al. ( 2006), while the dotted line and the dashed line indicate the 
observed range. We show three different characterizations of the host mass, as discussed in the text: (left) M sg the mass of both stars and gas; (middle) M lot the 
total mass; and (righ t) Mggidisk the total m ass within the half-mass radius. The last definition produces a scaling relation close to that observed between SMBHs 
and stellar spheroids ( Ferrarese et al. 2006), suggesting that it may be the closest analog of the mass M vi \ used by Cappellari et al. ( 2006) to characterize bulgeless 
spirals. Note that the two LT models in thin ci rcles , SG160-4 and SD220-4 do not fully satisfy the Jeans criterion, but judging from our resolution study, likely 
represent upper limits to the CMO mass. See j) l3.2l for details. 



particles is smaller than the Jeans mass) due to computation 
limit. However, from the resolution study in Figure Q] the 
CMO mass converges to roughly a factor of two even in runs 
where the gas particle mass is smaller than the initial CMO 
mass. We believe that these models represent the upper limits 
of the CMO mass of these two models. 

Figure [2] shows that for any definition of the galaxy mass, 
the mass of CMO correlates strongly with the mass of its host, 
despite the lack in our models of any explicit feedback from 
the CMO. It also appears that the correlation has a similar nor- 
malization and linear slope to that found for both star clusters 
and SMBHs. 

In reality, the neutral gas velocity dispersion in individual 
galaxies typically has a range of 6-15 km s" 1 so the scat- 
ter in the Mavio-A^gai relation, represented by the combina- 
tion of the LT and HT runs, may reflect the range of effective 
temperature we use here. Because the masses of the CMOs 
depend on the effective sound speed of the gas, it is desir- 
able to have a wide range of models with different effective 
sound speeds. Nevertheless, from the models we have sim- 
ulated, there is reasonable overlap of the two (LT and HT) 
temperature regimes, as shown in Figure[2] This suggests that 
the McMO~^gai relation is similar for different temperatures, 
apart from a factor of w 2 shift in the overall normalization, 
independent of host mass. 

It would be interesting to compare these findings directly 
with the observations, but su ch a comparison is difficult at 
present. IWalcher et al.l d2005l) have determined the masses of 



nuclear star clusters in 9 bulgeless spiral galaxies, and found 
them to range between 8 x 1 s M Q and 6 x 1 7 M Q . These ap- 
pear a factor of several smaller than the masses of the CMOs 
we find in galaxies withM tot = 10 11 - 10 12 M Q , but the masses 
of the hosts in their sample are not discussed, and their sample 
may also suffer from selection biases. 

Other comparisons are possible, but have to be interpreted 
with caution, since it is not clear how to relate the bulgeless 
disks in these models to the spheroid components that are 
found to correlate with CMOs in other observations. Nev- 
erthel ess, we present here two such brief comp a risons : 

(1) iBalcells et all d2003l) and iRossa et all d2006l) study 
late-type spirals with bulges. They find that the masses 
of the CMO at fixed bulge luminosity is about 3 times 
above the mass expecte d from the Mbh vs. Lbuige relation 
dMarconi & Huntll2003l) . In the left panel of Figure |2] we 
find a similar offset (relative to the MBH-^buige relation), if 
the total baryon mass is considered as a proxy for the bulge 
mass. One possible rationale for this comparison is that most 
of the baryon content of our simulated bulgeless disk galaxy 
may represent the stars making up the bulge of any galaxy that 
subseq uently form s out of t his bulgeless spiral . 

(2) ICote et alJ d2006l) : iFerrarese etail (|2006) and 
IWehner & Harris! J2006) study early-type galaxies, and 
find a correlation between the masses of the CMOs and the 
host galaxies. In particular, Fe rrarese et al. define M ga i, 
following the formula given by Capp ellari et al.1 d2006l) . as 
A^gai = aR e <j 2 /G, where G is the gravitational constant, a 
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FIG. 4. — Stellar mass of the galaxy as a function of total galaxy mass 
Mtot- The top panel shows the mass of newly formed stars (30% of the sink 
particles formed in the simulation, aside from the CMO), while the bottom 
panel shows the mass of all stars, including both newly formed ones and the 
old disk stars. 



FIG. 3. — This figure shows how the CMOs gain mass, as a function of 
total galaxy mass M tot . The top panel shows the ratio of the accreted mass 
Mace to the initial collapsed seed mass Mq, while middle and bottom panels 
show Mace and mean accretion rate of the CMO, respectively, as functions of 
Mtot of the galaxy. The mean accretion rate is the average over 1 Gyr (or the 
maximum simulation time if shorter) when Mcmo is measured, as explained 
in the text. The symbols are the same as in the previous plots, while the solid 
line is the least-square fit to the data. 



is the observed velocity dispersion, and a = 5 is a constant. 
This typically represents ~ 0.25 of the total galaxy mass (see 
below). Our definition of M ga i di s k is designed to mimic this 
quantity, and, as the right panel of Figure [2] shows, we find, 
within errors, that our bulgeless disk galaxies obey the same 
correlation. 

Both of the above comparisons suggest that the McMO-M„ d \ 
relation established in our simulated bulgeless disk galaxies 
may evolve into the observed correlation between Mcmo and 
their spheroid host, once the bulgeless disk galaxy (or a frac- 
tion of its mass) evolves into a spheroid. Establishing a direct 
connection, however, is not possible without modeling in de- 
tail the further evolution of both the host galaxy and the CMO 
(the latter, of course, may also gain mass). We next turn to 



CMOs that form in our simulations of merging galaxies. 

3.3. Toward a Physical Explanation of the Mcmo - M ga i 
Correlation 

The physics responsible for the mass correlation, in general, 
must be tied to the formation of a seed CMO, and its subse- 
quent growth by accretion. To investigate this issue further, 
we define the initial collapsed seed CMO mass as Mo, and 
the accreted mass as M acc . The final CMO mass is the sum 
of these two, Mcmo = Mo + M acc , where Mo is essentially the 
local Jeans mass. Our sink particles form at a fixed, arbitrary 
density n s j n k, and our models have a sound speed chosen from 
one of two values, so unsurprisingly, CMOs in galaxies with 
the same sound speed have a similar value of Mo. However, 
we find that M acc ^> Mo; that is, subsequent accretion plays a 
much more important role in determining the final mass of the 
CMO than initial collapse, as shown in Figure[3](fo/? panel). 

Figure [3] shows that M acc and the average accretion rate of 
the CMO increase on average with the total mass of the galaxy 
M tot (middle and bott om panels). W e expect the CMOs to 
accrete at roughly the Bondi (1952) rate (eq. Q]), so a more 
massive or gas-rich galaxy that has higher gas density, which 
enables higher accretion rate, ought to have higher values of 
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FIG. 5. — Correlation between the mass of the CMO and the global star 
formation efficiency in the galaxy in different models. Symbols are the same 
as in FigurefJ] The black lines show a least-square fit to the data points from 
models of each temperature. Because the CMO mass depends on the effective 
sound speed, so the lower Jeans masses in the LT runs result in a higher global 
SFE, and a smaller CMO. Our choice of two discrete values of c s represents 
the range of scattering (see text for more discussions). 



McMO 

Indeed, at fixed CMO mass and sound speed, equation fl} 
shows Mbh oc p. In a rotating disk in hydrostatic equilibrium, 

the ch aracteristic density scales as p cx M^ 3 (IWood & Loebl 
2000), assuming that the disk mass and scale radius are fixed 

1 /3 

fractions of the total galaxy mass M tot and size R cx M to , . 

2/3 

Hence, we expect M BH cx M to ' t , close to the power-law in- 
dex of 0.6 seen in the bottom panel in Figure[3] On the other 
hand, the gas-depletion time-scale tj, which halts both star- 
formation and the growth of the CMO, depends weakly on 
Mg a i (see Fig. [TJ. This leads us to a qualitative understand- 
ing of the McMO _ A^gai correlation seen in Figure [2] Since the 
CMO mass is dominated by accretion, we can very roughly 

2/3 

approximate it as Mcmo — Mbhtj cx M tot , while M ga i,di s k is 
roughly dependent on the stellar mass of the galaxy. As shown 
in Figure H] M star increases linearly with M tot , M star cx M tot . 
This generally supports the observed positive correlation of 
McMO-Mgai- Although this argument does not fully explain 
the observed linear relationship, it does ignore several factors 
that likely depend at least weakly on M tot , including the dark 
matter contribution to the galactic disk. 

The final mass of the CMO correlates tightly with the global 
star formation efficiency of the galaxy, with a normalization 
dependent on the effective sound speed, as shown in Figure[5] 
This further demonstrates a close link between the growth of 
the CMO and the build-up of the host galaxy through star for- 
mation. The dependence of the CMO mass on the effective 
sound speed can be understood from Bondi-Hoyle accretion. 
The initial mass of the CMO is close to the Jeans mass of the 
collapsing gas, Mcmo cx Mj cx p~ l l 2 cl, so M BH cx c\. Because 
of this strong dependence on the sound speed, our choice of 
two discrete values of c s at extremes of the expected range 
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FIG. 6. — History of CMO (solid) and total collapsed (dotted) mass (top) 
and accretion rate (bottom) in all five galaxy mergers. Total mass collapsed 
counts mass in all sink particles, a quantity related to the total amount of star 
formation in the galaxy by the local star formation efficiency. 



produces the discrete behavior seen in Figure[5] In reality, we 
would expect a scatter between these extremes. 

4. CENTRAL MASSIVE OBJECTS IN GALAXY MERGERS 

Previous simulations have demonstrated that major 
mergers of spiral galaxies genera lly trigger starbursts and 
produce elliptical galaxies (e.g.. | Tp omre & Toomre 1 1972 
Hernquisd [19891: IBarnes & Hernquistlll992t lHernquistlll992 
Mihos & Hernq uisl 1 19961: Dubinski, Mihos, & Hernquist 
19991: ISpringel(l200d IBarnes! l2002t iNaab & Burked 120031: 
Li et alj 120041) . Recently, simulations of galaxy mergers 
that include SMBHs have been carried out to investigate 
the SMBH-spheroid relations. Although some dissipation- 
less simulations showed that the M^-a relation could be 
maintained for mergers with modest energy and angula r 
momentum dNipoti et al J 120031: iBovlan-Kolchin et alj 120051) . 
it was pointed out that gas dissipation and star form ation are 
essenti al to this correla t ion j^ azantzidis et al. 2005). In par- 
ticular, Di Matt eo et alJ ([2005) suggest that thermal feedback 
from an SMBH suppresses both black hole accretion and 
star formation, giving the M BH -er correlation in ellipticals. 
Together, these studies suggest that the SMBH-spheroid 
correlations are closely connected to star formation, and that 
it may play a dominant role in determining the masses of 
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FIG. 7. — Surface brightness profiles of all five merger remnants at t he end 
of the simulations indicated in Figure[6] All five models fit well to a (Sersic 
I1963l[l9 68) profile with m = 4 (black line), as explained in the text. The 
legend gives the effective radius R e of each model. 




FIG. 8. — The relations between CMO mass and host mass and velocity dis- 
persion at the conclusion of our galaxy merger models. The value of Mcmo is 
the final mass of the CMO after completion of mergers, the a is taken at the 
mean values, and the error bar indicates the range over 10 4 viewing angles. 
The calculation of M ga i is described in the te xt. The black curve in the left 
panel indicates the best fit to observations by Tremaine et al. (2002), while 
those in the right panel are from observations by Ferrarese et al] 1200611 . 



SMBH and spheroids. 

In previous simulations, SMBHs are initialized in the pro- 
genitors following the M Bli -a correlation. In our merger sim- 
ulations, CMOs form dynamically as the galaxies evolve, in 
the same manner as in our models of isolated spirals. The 
growth histories of the CMOs in our merger models are shown 



in Figure [6] Because only the most massive CMO is consid- 
ered in the merger remnants, the CMO in galaxy mergers is 
identified as the most massive one. Similar to the isolated 
disk galaxies, the mass curves of the CMO increase rapidly 
in the initial phase. However, the detailed accretion history is 
modulated by the interaction between the progenitors. Close 
encounters trigger episodes of very high accretion rate, both 
on local and global scale, as shown in the bottom panel of 
Figure[6] The mass curves saturate after the completion of the 
mergers. 

In order to investigate the Mcmo-c correlation between the 
CMOs and the merger remnant velo city dispersions, we made 
simulated observations, following Gebh ardt et al.l (l2000t) . We 
measured the surface brightness-weighted, line-of-sight stel- 
lar velocity dispersion a of the spheroid within an aperture of 
size i? s /8, where R? is the effective rad ius. As mentioned in 
iBovlan-Kolchin, Ma. & Ouataertl d2005l) and lRobertson etalj 
(2006a), a standard technique for obtaining R e is to fit the 
projected surface b rightness profile, I{R) = L(R) / (4irR 2 ), to a 
(ISersidfT963lll968l) profile and derive the half-light radius: 



I(R) = I(R e )exp{-b(m)[(R/R e )^ m - 
1 /3 + 4/(405m), 



1]} 



as 



(2) 

given by 



where b(m) » 2m 
ICiotti & Berthl (1 19991) 

Bulges and early-type galaxies have surface brightness pro 
files I(R) that are usually well-app roximated by the de Vau 
couleurs law (Ide Vaucouleursl[l948l) 



logI(R)/I(R e )<x-(R/R e ) 



1/m 



(3) 



with Sersic index m = 4. Figure [7] shows the surface bright- 
ness profiles I(R) of all three of our merger remnants, assum- 
ing a fixed mass-to-light ratio of three. Th ese profiles ap- 
pear t o follow the m = 4 de Vau couleurs law (|de Vaucouleursl 
1 19481) closely, as also found by iHernquisl (1 19921) who stud- 
ied mergers of bulg eless spiral progeni t ors. W e note that 
iBalcells et aT] d2003l) and iGrahametaTI d2005l) proposed a 
smaller index, m ~ 1.7-3.5 for bulges, and it was suggested 
that only ellipticals with Mr < -21 mag have m ~ 4 profiles 
(e.g, [Graham & Guzman 2003). However, the derived value 
of m from fitting observations appears to depe nd on the data 
set and the observed wavelengths. As noted in Balcells et al. 
(l2003h . data from HST seems to give a much smaller m (e.g., 
m ~ 1.7) than that from Spitzer (e.g., m ~ 4-6) for the same 
sample, a nd low Sersic in dices are not expected from galaxy 
mergers. Graham (2006) reported a correlation between the 
masses of CMOs and th e Sersic index o f the h ost spheroids in 
the sample compiled bv lFerrarese et al.l d2006l) . More merger 
simulations with different initial conditions, mass ratios and 
orbital parameters are necessary to explore this relation. 

The effective radius R e is then derived from the half-light 
radius of each spheroid, as given in the legend of Figure [7] 
The mass-weighted line-of-sight stellar velocity dispersion 
a within R t j% is then calculated for each remnant galaxy 
for 10 4 random viewing angles. Figure [8] shows the result- 
ing correlation of Mcmo _(T f° r tne merger remnants. It ap- 
pears to agree very well with the obser vational correlation, 
M B h = 1-5 x 10V/200km s- 1 ) 4 - 02 M w dFerrarese & Merrittl 
2000: lTremaine et al.ll2002h . 

Figure [8] also shows the correlation McMO _ A^gai between 
the mass of CMO and M ga i of the host for the galaxy merg- 
ers. Here M ga i is calculated using the virial mass, follow- 
ing [Cappellari et al. (2006), M ga i = aR e a 2 /G, where G is 
the gravitational constant and a = 5. This is the best-fitting 
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virial relation based on the observables R e and a. The simu- 
lation agre es very well with the observed correlation Mcmo ~ 
0.002M ga i dFerrarese et alj2006l) . within the 1 a observational 
uncertainty, as indicated by the dotted lines in this figure. 

5. DISCUSSION AND SUMMARY 

Our simulations of gravitational collapse of gas in galax- 
ies enable us to model star formation in galaxies and follow 
the growth and evolution of CMOs. Our approach d iffers 
from previous studies such as Kazantz fdis et al.l (120051) and 
iDi Matteo et al.l (|2005) in the following ways. 

1 . Our resolution is sufficient to fully resolve gravitational 
collapse. 

2. Absorbing sink particles are used to directly follow and 
measure the mass of gravitationally collapsing gas. 

3. In our model, the black holes, or CMOs, are not set up 
a priori according to the M^-a correlation as in pre- 
vious approaches. Instead they form dynamically from 
gravitational collapse of gas, the same as star clusters, 
with both being represented by sink particles. In these 
simulations, the CMO is just the most massive sink par- 
ticle, which always is found in the galaxy center at the 
end of our simulations. 

4. Our model does not explicitly include feedback from 
either star formation or black hole formation. However, 
feedback from star formation is implicitly represented 
by our isothermal equation of state with relatively high 
effective sound speed of 6-15 km s" 1 . 

We emphasize that our models do not include explicit feed- 
back, magnetic fields, or gas recycling. However, we believe 
each will have minor effects on the final mass of the CMO. 
The assumption of an isothermal equation of state for the gas 
implies substantial feedback to maintain the effective temper- 
ature of the gas against radiative cooling and turbulent dissi- 
pation. Real interstellar gas has a wide range of temperatures 
but the root-mean-square (rms) velocity dispersion generally 
falls within the range 6-15 km s (e.g.. lElmegreen & Scald 
120041: iDib etaT1l2006l) . 

At least three effects may conspire to maintain the veloc- 
ity dispersion in the observed narrow range, as reviewed in 
iLi et al.l (1200 5b). First, radiative cooling drops precipitously 
in gas with sound speed < 10 km s" 1 as it becomes increas- 
ingly difficult to excite the Lyman a line of hydrogen. Sec- 
ond, direct feedback from the starburst may play onl y a minor 
role in quench i ng su bsequent star formation (e.g. iKravtsovl 
120031: iMonacol 120041) . perhaps because most energy is de- 
posited not in the disk but above it as superbubbles blow 
out (e.g. lFuiita et al.ll2003Ude Avillez & Breitschwerdtl2004l) . 
The energy input E from supernovae at the observed Galac- 
tic r ate drives a flow with rms velocity dispersion of ~ 9.5 km 
s -1 dde~ Avillez 2000), and the rms velocity dispersion depends 
only on E 1 / 3 (iMac Low! [l999t iMac Low & Klessenl 120041) . 
so a wide range of star formation rates leads to a narrow 
range of velocity dispersions. Third, magnetorotational in- 
stabilities may maintain a modest velocity dispersion even in 
the absence of any other feedback (Sellwood & Balbusl lT999t 
iDziourkevitch. Elstner. & Rudigeril2004l) . 

Radiative feedback from the central SMBH in massive 
galaxies will only Compton heat the gas much above 10 4 K if 



it the gas is sufficiently ionized. ISazonov et al.l (J2005) shows 
that this will occur if the ionization parameter 



= 1.4 x 10 s 



Mr 



1 cm 



10 8 M m n 



1 pc 
r 



(4) 



exceeds £ cr i t ~ 10 3 , with the Compton temperature of ~ 10 7 K 
being reached only for £ ~ 10 5 . In this equation, Mrh is the 
SMBH mass, L is the bolometric luminosity, which is as- 
sumed to be a constant fraction one-tenth of the Eddington 
luminosity, n the number density, and r the distance from the 
SMBH. If we examine the conditions at the Bondi-Hoyle ac- 
cretion radius for the SMBH, 



Rrh — 



GM B 



: (4.5 kpc) 



Mr 



10 8 M f . 



lOkms" 



(5) 



we find that gas with a density of n > (9 x 
IP" 3 cm- 3 )(M RH /10 8 M w ) will remain at 10 4 K (see also 
ISazonov et al.l 12005). This may contribute to maintaining 
the sound speed at roughly 10 km s" 1 in dense gas further in 
that might otherwise cool. On the other hand, the heating in 
the immediate neighborhood of th e SMBH may still lead to 
unsteady accretion as modeled bv lCiotti & Ostrikei] d2001l) . 

We should also emphasize that gas recycling from dying 
stars is not include d in our simulations. As computed explic- 
itly, for example, in Hai man et al.l d2004t) . stars can return over 
40% of their mass in gas. However, since the recycled gas 
might be heated up or blown out by the feedback, it may not 
be accreted by the CMOs. Our models with different effec- 
tive sound speeds satisfy similar correlations between Mcmo 
and host galaxy mass (though not the same relation between 
CMO mass and star formation efficiency). Furthermore, as 
discussed in § 2.2, the mass loss due to galactic winds is only 
a small fraction of the initial gas mass, suggesting that the 
effect of feedback on the mass of the CMO may be minor. 

The simple model of sink particles that grow with ap- 
proximately spherical Bondi-Hoyle accretion reproduces the 
observed Mcmo-ct and McMO _ M ga i correlations reasonably. 
Both the accreted mass of CMO and the mass turned into 
stars scale with the total mass of a galaxy, offering a plau- 
sibl e explanation for the CMO-host relation. We should note 
that|Escala (2006) also reproduces the Mrh-c correlation by 
modeling accretion o f BHs from a Shakura-Sunyaev thin disk 
dShakura & Sunvaevlll973l) . This suggests that the accretion 
rates through these two different mechanisms may be sim- 
ilar, and that the final mass of the black hole or the CMO 
may be largely determined by the depletion of accretable gas. 
The fact that various previous approaches, including fee dback 
models ("e.g. IDi Matteo et al.ll2 005: Saz onov etaT]l2005l) . star 
formation model s (e.g.. [Kazaptzidis et alJ 120051) . and accre- 
tion models (e.g., Escala 2006), have more or less succeeded 
in reproducing the SMBH-host correlation suggests that gas 
depletion in the central region is crucial to this correlation. 
Our results support the idea that star formation plays the major 
role in gas consumption and determination of the final masses 
of the CMO and the stellar component. 

The close correlation that we find between the CMO mass 
and global star formation efficiency in the host galaxy sug- 
gests that the CMO-host link may be universal, in that it is the 
result of the co-eval growth and evolution of the CMO and the 
host galaxy, and therefore it does not strongly depend on the 
morphology, type, or mass of the galaxy. A systematic search 
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for CMOs in the nuclei of bulgeless disk galaxies would offer 
a test of this conclusion. 

In summary, our simple model of accretion and star forma- 
tion reproduces quantitatively the observed Mbh-c correla- 
tion for galaxy mergers, and suggests a universal McMO-^ai 
relation over a wide range of galaxy mass and morphological 
types. We find that the CMO builds up its mass through accre- 
tion, and that there is a direct correlation between the global 
star formation efficiency of a galaxy and its CMO mass. Our 
results suggest that star formation may play an important role 
in producing the fundamental mass correlation between the 
central massive objects and their host galaxies. 
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